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Abstract 

We present an analitical study of the dynamical process of the approach to steady 
state for a driven diffusive system represented by the microemulsion phase of a 
ternary mixture. The external applied field is given by a "plane Couette" shear flow 
and the problem is studied within the framework of a time-dependent Ginzburg- 
Landau model. A Fokker-Planck equation for the probability distribution of the 
concentration fluctuations is developed in a self-consistent one-loop approximation 
and solved exactly, giving an analytical expression for the dynamic correlation func- 
tions of the system. For comparison to experimental work, we also show grey-scale 
plots of the scattering function at different times during the dynamical process for 
two different shear rates. 
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1 Introduction 



Ternary mixtures made up of oil, water and surfactant agent are of consider- 
able interest for both fundamental issues as well as for applications [1-3]. In 
these systems, three-phase coexistence is commonly observed between a water- 
rich phase, an oil-rich phase and a microemulsion. The latter is a mesoscopic 
structured phase where coherent domains of oil and water on scales between 
100 and 1000 A form disordered isotropic intertwined structures on larger 
scales [4] . A very similar mesoscopic structure is observed in sponge phases in 
acqueous amphiphilic solutions [5,6] and in homopolymer blends containing 
small amounts of diblock copolymers which act as a surfactant [7-11]. 
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After the earlier work by Gompper and Shick [12] in which they demon- 
strate that the failure of the microemulsion to wet the oil-water interface 
is a consequence of the microemulsion internal structure, various theoretical 
advances on similar ternary systems have been made [13-16]; in particular 
Kielhorn and Muthukumar [16] have investigated the weak segregation be- 
haviour of a ternary mixture of A homopolymer, B homopolymer and AB 
diblock copolymer, showing several features which are identical to those seen 
in traditional oil/water/surfactant systems. These amphiphilic systems are 
now receiving renewed attention, since several light and neutron scattering 
experiments revealed the hallmarks of glassy dynamics, in the sense of large 
relaxation times, on their mesoscopic structures [17-22]. Other ways of es- 
tablishing non-equilibrium states with slow dynamics are also possible. For 
instance, microemulsions can be driven to a slow out of equilibrium station- 
ary regime by a weak shear. This perturbation accelerates the dynamics in 
such a way that the structural relaxation time decreases with increasing shear 
strength (a phenomenon usually known as shear-thinning). The shear flow can 
therefore be used as a probe of the nonequilibrium properties of these systems 
by adopting the shear rate 7 as a control parameter [23]. Several scattering 
experiments have shown the very complex and diverse effects of shear flow on 
these systems after the steady-state has been reached [24-32]: nonlinear effects 
due to the flow give rise to strongly anisotropic suppression of the composi- 
tion fluctuations; the phase diagram is affected by the flow which can induce 
order-disorder transitions; the rheological behaviour is highly non-Newtonian 
with significant intriguing shear-thinning and shear-thickening behaviour. 
In the following, we shall find useful to refer to the concentration fluctuation 
(equal time) correlation function, calling it "structure factor" , although such 
a name usually indicates the density fluctuation. 

Theoretical efforts to understand these complex problems have also been made 
intensively. A mean field theory able to describe the scattering function of a 
diblock copolymer in the presence of a steady homogeneous flow field was de- 
veloped by Fredrickson [33] as a generalization of the equilibrium theory of 
Leibler [34]. The theory predicts the wave vector dependence of the scatter- 
ing intensity and the distribution function of composition fluctuations in the 
disordered phase. The concept of an increase of the order-disorder transition 
temperature due to the shear was first studied for small shear rates by Cates 
and Milner [35] who applied a self-consistent Hartree approximation to the 
nonequilibrated system and solved for the shear effect on this Hartree term. A 
more general self-consistent approach, not limited to small shear rates and with 
predictions which are in good agreement with experimental results [36-38] , was 
developed by Huang and Muthukumar [39] who studied numerically the effect 
of shear rate on order- disorder and order-order transitions involving lamellar 
and cylindrical morphologies as well as the dependence of steady-state scatter- 
ing function on shear. A one-loop self consistent Ginzburg-Landau model for 
describing the steady-state rheological behavior of oil/water/surfactant mix- 
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tures under a plane Couette shear flow was considered by Patzold and Dawson 
[40] who gave a numerical evaluation of the structure factor together with an 
analytical analysis in the limit of vanishing shear rate 7. Their model was 
exactly solved in [41] for any value of 7 showing that in a specific shear range 
the structure factor is characterized by four pronounced peaks in the plane of 
the shear and flow directions. This phenomenon, which was already predicted 
by Drolet et al. in a slightly different model [42] and also reported in scatter- 
ing experiments on segregating mixtures [30,25], was interpreted in terms of 
a complex spatial pattern where interfaces with two different orientations co- 
exist, with relative abundance tuned by 7. On the other hand, the relaxation 
dynamics of microemulsions towards a steady state, which can be extremely 
slow, is also very intriguing: in fact, the absence of time translation invariance 
gives rise to a complex behaviour where different modes relax on different time 
scales and differently compete with the hydrodynamics time scales. 

A systematic study of the dynamical behaviour of microemulsions and sponge 
phases in thermal equilibrium and in the presence of hydrodynamic modes 
has been developed by Hennes and Gompper in [43]: they use a field-theoretic 
perturbation theory which makes use of "response fields " [44-46] and whose 
key ingredient is a self-energy matrix having a block structure identical to 
the Gaussian correlation and response matrix. Within that framework, they 
were able to calculate the frequency-dependent (complex) viscosity, the sound 
velocity and damping in a one-loop approximation. However, the macroscopic 
relaxation of the system and the role of the different time scales that are 
involved are most clearly shown by the evolution of the scattering intensity 
and, in particular, its asymptotic behaviour for short and large times. This 
was calculated in [43] by using the first Born approximation for evaluating the 
correlation matrix and, as it was already pointed out by those authors, it is not 
clear to what extent such behaviour relies on the adopted approximation. As a 
result, a self-consistent calculation of the dynamic behaviour of microemulsions 
and sponge phases is expected to clarify this important point. 

Of concern in this paper is the microemulsion phase only and its slow relax- 
ation dynamics towards the steady state, associated with composition fluctu- 
ations under a plane Couette shear flow. We will focus on the time- dependent 
static structure factor (i.e. the average of the concentration fluctuation at equal 
time, whose time dependence comes from the transition from equilibrium to 
steady state), and perform a self-consistent one-loop calculation within the 
model previously adopted in [41], which has the significant feature that results 
can be formulated in terms of the physical correlation lengths, thereby making 
direct comparisons to experiment possible. The aim of the present work is to 
contribute in making a contact between the mesoscopic structural character- 
istic of microemulsions and the general macroscopic dynamic features of soft 
materials. The plan of this paper is as follows. In Sec. 2 we present the adopted 
Ginzburg-Landau model for the microemulsion and derive the Fokker-Planck 
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equation for the probability distribution of composition fluctuations. In Sec. 3 
we obtain its dynamic solution which gives an analytical expression for the 
time-dependent static structure factor. We discuss the scattering function con- 
vergence towards its steady-state value and show its morphologies during the 
relaxation process. In Sec. 4 we draw our conclusions. 



2 Ginzburg-Landau Model 

We shall consider a binary fluid of oil and water at 50:50 concentration, with 
the addition of non-ionic surfactant agent, in thermal equilibrium at a temper- 
ature T > T c . In order to describe the ternary mixture, we adopt a Ginzburg- 
Landau model with a single scalar order parameter <p representing the lo- 
cal concentration difference between oil and water components. The following 
coarse-grained Ginzburg-Landau free-energy, firstly proposed by Gompper and 
Schick [3], has been used extensively for investigations of the static structure 
and phase behavior of these systems [47,3,12]: 



where /(</?) is a polynomial function of if which has three minima <fi < (fi2 < 
</? 3 , with ip2 — 0, corresponding to the three homogeneous bulk phases, oil- 
rich, microemulsion and water-rich respectively. Moreover g{f) = go + g2 ( p 2 , 
with g < and g 2 > 0, in order to grasp the experimental evidence that the 
peak in the structure factor in the oil and water phases is at zero wavevector 
(so that g((p) must be sufficiently positive there), while the peak is at non- 
zero wavevector in the microemulsion (so that g{f) must be negative there). 
The last term, with a coefficient c > 0, makes stable the free-energy at large 
momenta and weights the curvature of the interfaces where surfactant gathers. 
As we are concerned with the microemulsion phase, i.e. the disordered but 
structured region of the phase diagram where a > and a stability condition 
Aac — gl > is fulfilled, we will use a piecewise form of /(</?): 



as proposed by Patzold and Dawson [40], where the quartic term has been 
included in order to capture the fluctuation effects in the critical region. Al- 
though local concentration of surfactant does not appear in the free-energy, 
its average value is taken into account through the coefficients go and c: in 
particular the negative value of g , known as the surface tension parameter, 
allows us to model a characteristic feature of the microemulsion phase, i.e. 






(2) 
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the tendency of surfactant to create spontaneously interfaces. More refined 
models should include an additional order parameter to describe the density 
of surfactant (for a review about this point see [3]). In the following we will 
be concerned with the response of the system to a plane Couette shear flow 
where the average velocity field is expressed as 

v(x) = A /ye x (3) 



where 7 is the shear rate and has the dimensions of a frequency. The reference 
frame adopted is such that the x-axis is along the flux, the y-axis is along the 
velocity gradient and the z-axis is along the vorticity direction. e x appearing 
in Eq. (3) is the unit vector in the flux direction. We will only consider shear 
flows such that the molecular relaxation times are less than the inverse of the 
shear rate 7, i.e. the Deborah number D = « 1. Moreover, we will assume 
that, due to the highly viscous nature of the mixture, the fluctuations of the 
velocity field are small and that their coupling with the fluctuations of the 
order parameter is negligible [33,48]. Under these assumptions it is legitimate 
to apply the stochastic model given by Onuki and Kawasaki [49] and study in 
the momentum space the composition fluctuations by the following Langevin- 
like equation 



where the Onsager coefficient T is assumed to be a constant at a temperature T 
sufficiently higher than T c (this assumtion cannot be extended to homopolymer 
mixtures where a k-dependent Onsager coefficient must be used [50]) and 
SF/Sip is the usual thermodynamic force. 77 is a stocastic white noise describing 
thermal fluctuations into the system: it is Gaussian distributed with mean 
zero and variance constructed to satisfy a (generalized) fluctuation-dissipation 
theorem [51]: 



MM)) = o (5) 

(r]{k, t) rj(k>, t')) = 2 T k 2 5{k + k>) 5{t - t') (6) 

As usual, it is assumed that the fluctuating force is uncorrelated in time, 
which implies that one considers times much larger than the characteristic 
times during which internal modes equilibrate. Thereby, the stochastic process 
becomes Gaussian and Markovian [52]. In order to linearize the problem, we 
further assume that it is legitimate to approximate the term / 8tp(—k, t) 

as x(k) vikyt), where x 1 (^) can be regarded as a static susceptibility that 
depends on the correlation functions evaluated in the steady state. A self- 
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consistent one-loop approximation gives 

x (k) =a R + g R k 2 + ck* (7) 



dk 

S? = p* I -—\k\ n C{k,t) (8) 



with a R = a + 6/2 + g 2 S%°, g R = g + g 2 S%°, and 

dk 

|fe|<A 

In Eq. (8), C(k,t) = (<p(k,t)<p(—k,t)) is the time-dependent static structure 
factor and A is a cutoff proportional to the inverse of the length of the surfac- 
tant molecule. It is worth to note that, apart from the presence of renormalized 
coefficients, x~ l {k) preserves the same functional expression of the equilibrium 
scattering function firstly deduced by Teubner and Strey in their seminal pa- 
per [47]. Such a form yields for a > 0, g < and c > a single broad scat- 
tering peak at a non-zero wavevector, and a k~ A decay at large k which are 
properties experimentally observed for a variety of microemulsions containing 
comparable amounts of water and oil [53-56]. The present Hartree approx- 
imation has been used by Huang and Muthukumar (see [39] and references 
therein) for studying the steady-state structure factor in block copolymers. 
In the nonequilibrium problem in which we are interested in this paper, it 
will be used to study the dynamical process of the approach to the steady 
state. It represents the main simplification in our model and it is deemed to 
be adequate only close to the steady state (see also [57] and references therein, 
where a similar argument was used to study the approach to the steady state 
of staggered non-equilibrium particle systems and Ising chains). Therefore, the 
linearized Langevin equation for the compositions fluctuations is written as 

^ <p(h, t)=^k x -^-<p(k, t)-Tk 2 x (k) <p(h, t) + n(%, t) (9) 
The set of Eqs. (5)- (9) is equivalent to the following Fokker-Planck equation 

— * — * 

for the probability distribution of the (p(k), P((p(k),t) [58]: 



+ r\k\ 2 (x(k) <p(k) + P(v(k),t) (10) 

Eq. (10) completely defines the stochastic relaxation process in which we 
are interested and represents the starting point for the analysis of the time- 
dependent static structure factor that will be presented in the next session. 



6 



3 Results and Discussion 



Eq. (9) states a linear dependence of ip(k, t) from the noise r)(k, t): as the noise 

— * 

is Gaussian distributed, then the probability distribution P[ip(k, t)] must be 
Gaussian at each time. In fact the steady solution of Eq. (10), which represents 
the probability distribution for the field ip in the steady state, is given by 



JetC (k) 




1 r dk 

2 J 72tt 



dk (p(k) <p(—k) 



2 J (2tt)3 Co (k) 



where Co stays for the steady-state structure factor and has the following 
expression [41] 

oo 

Co(jfe) = Jdt e~fi " {m)dC T | k(t) | 2 (12) 



where 



kit) = (k x , ky + ^1 , kz ) (13) 



and 

uj(k) = Tk 2 X (k) (14) 

is the thermal decaying rate, i.e. the rate at which order parameter fluctuations 
having a wavelength A = 2n jk relax at thermal equilibrium (in the absence of 
flow) [59,49]. Eqs. (7), (8) and (12) form a self-consistent relationship which 
accounts for the fluctuations under shear. By consideration of Eq. (12) it is 
apparent that for 7 >> uj{k) the composition fluctuations will be strongly 

— * 

distorted by the flow, i.e. the flux forces the C (k) to be highly anisotropic. 
The integral in Eq. (12) is clearly positive, as it should be for a variance, and 
it is possible to verify that the solution (12) is unique: in fact every other 
solution different from that does not recover the Ornstein-Zernike form in the 
limit 7^0 [49]. It is worth to note that C (k) does not depend on the initial 
conditions. It will be shown in the following that the steady-state structure 
factor can be obtained from the time-dependent static one through a limit 
calculation (t — > 00) after which the correlation loses the memory of its initial 
value. Because of the Markovian character of the stochastic process and the 
translational invariance (see [49]), we make the following ansatz for the time- 
dependent probability distribution of the field (p during the dynamical process 
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of the approach to steady state, that will be completely proved a posteriori: 



P(<p(k),t) = G(t)exp 



dk (fi(k) <f(—k) 



(2tt) 3 C(k,t) 



(15) 



where G(t) = [l/ (det C(k,t)) 



ll/2 . 



is a normalization factor, and C(k, t) is the 
time-dependent static structure factor, which is given by the time-dependent 
solution of the following equation 



8c(M)_. t a^_ M , )cfti) + 2 „ 2 



dt 



dk,, 



(16) 



Eq. (16) can be obtained formally by taking the time derivative of C(k,t) = 
(ip(k,t) ip(—k,t)) and using the linearized Langevin Eq. (9), provided a regu- 
larization procedure is used for dealing with the difficulties coming from the 
fact that time differentiation and averaging do not commute. Such a regular- 
ization entails defining (ip(k, t) i](—k, t)) = \im t i^ t - (ip(k, t) r](—k, t')) 6(t — t'), 
where the presence of the step function 0(t — t') reflects the causality of the 
Langevin equation, and then setting this equal to | (r)(k,t) r)(—k,t)) [60]. 
The first and second term in the R.H.S. of Eq. (16) are the convection term 
and the diffusion term respectively, while the third term is the fluctuation 
source. When the shear rate is large enough, the hydrodynamic force strongly 
interferes with the thermal decaying process and a nonlinear effect of suppres- 
sion of the composition fluctuations takes place. We shall begin by solving the 
homogeneous problem: 



ec(k,t)_.^ecg£ = C{U 



dt 



dk y 

C(k, 0)=C mit (k 2 



(17) 
(18) 



where C ini t(k 2 ) = 1/ (a + g k 2 + ck A ) is the well known isotropic scattering 
function of the microemulsion at initial time, in the absence of flow [47,3], and 
its projection on each Cartesian plane is that of a circular volcano with radius 
proportional to the inverse correlation length of the system, as it is shown in 
Fig. (1). 

The solution of the homogeneous problem is found by using the method of 
characteristics and it has the following form: 



C hom (k,t) = C init (k 2 (t))exp 



- / uj(k{t-t'))dt' 



(19) 
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Fig. 1. The structure factor at initial time. 

where u(k) has been defined in Eq. (14). The full solution of Eq. (16) with 
the initial condition given in Eq. (18) is 

t 

C(k- 1) = C hom (k- t)+fr\ k{m) | 2 e--T a,( * (c))dC tim (20) 

o 



Taking the limit t — > <x>, it is straightforward to verify that C(k,t) — > Co(k), 
i.e. the time-dependent field distribution P(ip(k),t) converges (in a weak 
sense) at infinite time towards the steady-state distribution P (<p(k)): this 
justifies a posteriori the previous ansatz (see Eq. (15)). Eq. (20) is the prin- 
cipal result of this paper: it represents an analytical expression for the time- 
dependent static structure factor of microemulsions in the presence of a plane 
Couette shear flow. We shall now consider the behaviour of the structure factor 
for short times. In this limit we have 



c{k-t) = c hom (k- o) + r k 2 i - c imt (k 2 (t)) x (k) 



t + 



+ dC ™^ 2) jk x k y t + 0(t 2 ) (21) 

i.e. the scattering intensity grows linearly in t. A linear dependence on time t 
has been originally predicted by Hennes and Gompper in [43] within a differ- 
ent model. Nevertheless, the present self-consistent calculation differs from the 
former one in two important respects: firstly, the expression given in Eq. (21) 
shows a dependence on mode-coupling terms through x(fc) that does not ap- 
pear in [43]: their presence is likely to be important in predicting the cor- 
rect approach to the steady state. Secondly, the third term on the R.H.S. 
of Eq. (21) explicitly breaks the isotropy of the short time behaviour of the 
scattering function, unless k x k y = 0. The significance of this anisotropic 
term is apparent in the following analysis and is amenable to experimental 
test. The predictions of our model are now investigated on the strongly struc- 
tured microemulsion whose steady-state characteristic features have already 
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been studied in [41] (to which we also refer for details on the specific values 
adopted for the system parameters). The role of the shear rate 7 as a control 
parameter for probing the nonequilibrium properties of microemulsions has 
been already discussed in the Introduction. Therefore, we will now use two 
distinct values of the shear rate (7 = 2 and 7 = 8) which, even if both belong- 
ing to the region in which shear-thinning behaviour is observed, still produce 
very different mesoscopic structures into the system. The characteristic mor- 
phologies of the time-dependent static structure factor projections on both 
the planes k z = and k y = 0, are shown in Fig. (2) (for 7 = 2) and in Fig. (3) 
(for 7 = 8) at four different times. The first time (t = 100) is chosen to 




Fig. 2. Projections of the time-dependent static structure factor on the planes 
k z = (left column) and k y = (right column) at different times t=100, 300, 500, 
1000 (from top to bottom). The value of the shear rate is 7 = 2. 

be sufficiently small to observe the early deformations of the initial isotropic 
volcano shape, while the last one (t = 1000) is such that the steady-state 
morphologies previously found in [41] are recovered. The effect of the flow on 
the shape of the structure factor is to stretch it along a line oriented at 135° 
to the fc x -axis and rotate it towards the k y -a,xis (fc^-axis) in the k z = plane 
(k y = plane). At any fixed time during the relaxation dynamics, the degree 
to which the scattering function is stretched and rotated depends on the value 
of 7 and grows with it. It is worth to stress that the short-time morphologies 
in the k y = plane are less anisotropic than those in the k z = plane as 
a result of the absence of the term oc 7 k x k y t in Eq. (21) when k y = 0: the 
predicted different degree of isotropy should be observable in scattering exper- 
iments and its possible confirmation would represent an important test of the 
present theory. Finally, the steady-state morphology of the structure factor 
strongly depends on 7: in the case 7 = 2 the predicted four peaks structure 
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Fig. 3. Projections of the time-dependent static structure factor on the planes 
k z = (left column) and k y = (right column) at different times t=100, 300, 500, 
1000 (from top to bottom). The value of the shear rate is 7 = 8. 

can be related to a biphasic region with different allignment while above a 
certain shear rate the microemulsion is predicted to transform into a lamellar 
phase oriented along the flow, as it is observed for 7 = 8. 



4 Conclusions 



In the present paper we have studied the response of microemulsion to the sud- 
den application of a steady plane Couette shear flow. Analytical expressions for 
the dynamic scattering function and the distribution function of composition 
fluctuations have been developed and their convergence to steady-state values 
proven. The prediction of the model have been investigated on a structured 
microemulsion whose time-dependent static structure factor has been evalu- 
ated and explicitly shown at different times during the relaxation dynamics 
towards the steady state. The present results can be tested in neutron or X- 
ray scattering experiments through which vital informations can be gained to 
further developing the analysis of the dynamical process. 
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